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An extension of tlie MSSM called the /iz^SSM does not allow a conventional thermal 
leptogenesis scenario because of the low scale seesaw that it utilizes. Hence, we 
investigate the possibility of electroweak baryogenesis. Specifically, we identify a 
parameter region for which the electroweak phase transition is sufficiently strongly 
first order to realize electroweak baryogenesis. In addition to transitions that are 
similar to those in the NMSSM, we find a novel class of phase transitions in which 
there is a rotation in the singlet vector space. 



1. INTRODUCTION 



An extension of the MSSM called the /ii/SSM [T] is a model similar to the NMSSM [2] 
(with the usual Z3 charge assignment) except that the singlet whose vacuum expectation 
value (VEV) gives rise to the /z term also serves the role of a right-handed neutrino, thereby 
violating R-parity. Because the VEV generates the /z-term and the right handed neutrino 
mass, the right-handed neutrino masses are of order TeV, leading to a low scale type I seesaw. 
Given the absence of a high scale seesaw, thermal leptogenesis is difficult in the /xi/SSM. 
Hence, it is interesting to consider whether or not electroweak baryogenesis (EWBG) [3] can 
occur in this class of models. One of the most stringent constraints of EWBG on the /iz/SSM 
is the requirement of a sufficiently strongly first order phase transition (SFOPT) such that 
the created baryons are not washed out [1]. 

Because the /xz/SSM contains 3 singlet chiral superfields (right handed neutrinos), mainly 
motivated by generality, standard model generation replication pattern, and phenomenologi- 
cal convenience [HEI, there is a "larger" SFOPT parameter space for EWBG when compared 
to the NMSSM. More precisely, there can be SFOPT where the singlet VEVs rotate in the 
singlet vector space during the electroweak phase transition. The price paid for this is a 
more complicated global minimum analysis at both finite and zero temperatures. The aim 
of this paper is not to uncover the most general parameter space consistent with EWBG, 
but is to simply give a couple of parametric regions to show the existence of possibilities. 

Depending on the path of the phase transition, the exact /xz/SSM parametric dependence 
of the phase transition strength v{Tc)/Tc is complicated. Nonetheless, we find that it is 
typically true that to achieve SFOPT, the parameters are close to satisfying the following 



* Electronic address: 
'''Electronic address: 



danielchung( 


§wisc.edu 


aj long® wise 


edu 





2 



condition: 

(1) 

where Ees is the effective cubic couphng, Aeg is the effective quartic couphng, and f (0) 
is the magnitude of the scalar field space VEV (including both the Higgs and singlets) at 
zero temperature. Physically, this corresponds to the parametric region where the critical 
temperature is small compared to w(0) during the electroweak phase transition. In the 
examples provided in this paper, whether or not the SFOPT proceeds from the origin, the 
leading nonvanishing value of Ees in the /xz/SSM arises from the soft terms 



E 



1 ^ 



(2) 



where i>f are singlet fields. The dimensionful coupling ax is distinguished from in that 
a\ also enhances the mixing between the Higgs sector and the singlet sector. The leading 
contribution to Acg comes from the superpotential and D-terms. 

Beyond these general results, we find a somewhat interesting feature because we focus on 
the parametric region analyzed by [5]. In this parametric region, an approximate S3 symme- 
try (permutation symmetry) arises due to the right handed neutrino generation independence 
of the non- Yukawa couplings and the smallness of the neutrino Yukawa couplings. Hence, 
to avoid any extra complications associated with domain wall formations, one might naively 
try to avoid §3 symmetry breaking phase transitions by considering parameters which yield 
zero temperature vacua preserving §3. Hence, this is the boundary condition that we impose 
in this paper. Interestingly, we find that despite this boundary condition, §3 is typically 
spontaneously broken multiply at finite temperatures in a way that is sensitive to quantum 
radiative corrections. As the temperature is lowered from high temperatures, this leads to 
multistep phase transitions starting from the trivially S3 symmetric vacuum in which all 
VEVs vanish. The electroweak symmetry breaking phase transition occurs with S3 symme- 
try restoration to a vacuum in which all sneutrino VEVs are identical and nonvanishing. 
We also find one step SFOPTs in which the scalar fields (including the singlet fields) make 
a transition from the origin to the electroweak symmetry breaking minimum. The numeri- 
cal values of the parametric regions uncovered in this paper is in the paragraph containing 



Eq. (46) and regions Ilia and Illb depicted in Fig. [5j 

There have been many studies of EWBG and the electroweak phase transition in models 
with gauge singlets; some of these are Refs. P-E], and we will discuss others throughout 
the remainder of this section. Since our work is most closely related to previous work on 
SFOPT in the NMSSM, we give here a little preview of some of the differences between our 
work and select previous works, in addition to the multidimensional aspect stressed above. 
In Ref. [H], SFOPT in the context of the NMSSM is first analyzed and the author points 
out that the tree level cubic term coming from the soft SUSY-breaking sector is important. 
Note that Ref. [T3] uses the definition of critical temperature in which scalar mass squared 
matrix develops a vanishing eigenvalue. We take a more robust definition of T^. being the 
temperature at which a new coexistence phase occurs even though this definition is harder 
to implement in practice. 

The authors of Ref. p,6j also analyze the NMSSM, but they include a /x-term on the 
basis that it is more general and its nonzero value eliminates the Z3 symmetry which can 
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be cosmologically dangerous with respect to the problem of domain wall formation fTf]- 
The nonzero /i-term leads to false vacuum not being at the origin. In this case the critical 
temperature criterion used by Ref. [15j is invalid. Therefore, the authors of Ref. [16i| take 
the coexistence phase definition of critical temperature as we do in this paper. They also 
include a bilinear soft term in the Higgs which breaks the Z3 symmetry. Although we do not 
include such Z3 breaking terms directly, we will assume that nonrenormalizable terms can 
be included to obtain acceptable phenomenology with respect to any possible domain wall 
formation. However, it is to be noted that Z3 breaking can often lead to UV instabilities in 
the singlet tadpoles, making the UV stability of these theories (including the one considered 
in this paper) a model building challenge as noted by [T7] . 

The analysis [18] considers the generalized NMSSM similar to [16]. They run 9 parameters 
with a popular choice of "universal" boundary conditions from the GUT scale down to the 
electroweak scale to generate their model. They do not reject metastable vacua based on 
the intuition that longevity of the false vacuum on the horizon scale today is not difficult to 
attain. To be conservative and to avoid potentially complicated discussions of met ast ability, 
we accept only stable vacua in this paper. 

A model related to the NMSSM and the /xi/SSM is the nMSSM in which the discrete 
charge assignment is modified as to eliminate the singlet cubic term in the superpotential. 
This model was analyzed by [19] for SFOPT. For a significant portion of the parameter space 
in which SFOPT occurs, a linear tadpole term in the superpotential plays a significant role 
in contrast to our scenario. 

The analysis of [20] considers the EWPT in an extension of the SM which adds a real 
singlet S. These authors find a large region of the parameter space of their model that is 
consistent with SFOPT and LEP Higgs search bounds. They argue that the strength of 
the phase transition can be enhanced by 1) choosing a large negative value for the SH^ 
coupling, 2) choosing a negative value for the S^H"^ coupling, or 3) allowing the singlet to 
have a nonzero VEV before the electroweak symmetry is broken. In the language of this 
paper, the first two points correspond to increasing E^s and decreasing Aes, respectively. 

Before we begin the main body of the work, let us list here all the caveats to our anal- 
ysis. We do not take into account explicitly the high energy Landau pole constraint (i.e. 
perturbativity up to the GUT scale) because we will take the attitude that the /xi^SSM is 
well motivated mainly by its ability to have all fields participate at low energy and thereby 
have potential measurability. Nonetheless, the parametric region that we uncover lies at the 
border of perturbativity up to the GUT scale (inferring from the work of Refs. [SJET]), which 
means that the UV cutoff for our theory can be taken to be far above the TeV scale. We do 
not take into account explicit Z3 breaking effects because a small amount of breaking can 
address most cosmological domain wall problems, as we later demonstrate. We do not take 
into account explicit CP violation effects in the phase transitions as this will typically lead 
to less than order 10% effects since CP violating phases compatible with phenomenology 
are typically order 0.1 or smaller. For robustness, we accept in this paper as phenomeno- 
logical possibility only absolutely stable global zero temperature vacua instead of analyzing 
the phenomenological possibilities of metastable vacua. Finally, all of our numerical work is 
kept in control to only order 10% accuracy. 

The order of presentation is as follows. In the next section, we present the Lagrangian 
including its discrete symmetry properties and radiative/thermal corrections. The section 
concludes by highlighting the yUz/SSM differences from the NMSSM scenario. In Sec. [3| 
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we describe the parametric region relevant for SFOPT in terms of one-dimensional field 
space slice parametrization. There we also qualitatively describe how the multidimensional 
paths of the phase transition and discrete symmetries play a role. Next, in Sec. |4| we 
explicitly show that singlets do not play a significant role in terms of numerical value of the 
sphaleron action controlling the B+L violating rate in the broken phase. The main numerical 
results are presented in Sec. [5] where the explicit existence of the SFOPT parameter region 
is demonstrated. Details of the transition paths organized in terms of discrete symmetries, 
phenomenological bounds placed, and explicit mass spectra for a sample parametric point are 
given. In Sec. |6} we demonstrate that the cosmological domain wall problem is easily evaded 
with an inclusion of a weak Z3 symmetry breaking operator in our scenario. We then conclude 
with a summary of the results. Several appendices then follow giving useful technical details. 
In Appendix A, we list the field-dependent mass matrices used for computing the effective 
potential. In the next appendix, we give details regarding the approximate thermal masses 
used in the paper. In Appendix C, we describe analytically the boundaries of in Fig. [Sjwhich 
is one of our main results. Finally, in Appendix D, we show that it is generically possible to 
construct a nonrenormalizable Z3 superpotential to obtain a CP conserving global minimum 
in the absence of any explicit CP violating parameters. 



2. THE THERMAL POTENTIAL DIFFERENCES BETWEEN THE NMSSM 

AND THE ;^zySSM 

The /iz/SSM that we consider in this paper is specified by the following superpotential 
and soft terms 

i ^ ^ 

(3) 

3 

i 

2 / 3 \ 3 

+ ^m^Ji/,p-M^M,AA + h.c. j + 

i \ i / i 

3 

+ Yl {auQi-H2Ui + adHi-Qidl + aeHi-licI + a^Li-H2vt + h.c.} . (4) 

i 

Where indicated by a dot, the SU(2) indices are contracted with the antisymmetric tensor 
and ei2 = 1. First, note in addition to the usual Z3 symmetry used to forbid an explicit jjL- 
term, there is an exact CP symmetry due to the reality of the coupling constants. We ignore 
the CKM phases since these will only give corrections smaller than the (9(10%) accuracy 
that we are aiming for in this paper. The CP transformation in the scalar effective potential 
effectively takes each scalar field to its complex conjugate. Next, note that the couplings 
of the sector to the SM were taken to be generation independent, except for the Yukawa 
couplings, and that the singlets do not couple with one another directly in the superpotential. 
This choice is motivated by trying partially to match the work of [5j. Hence, we see there is 
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an enhanced S3 symmetry (permutation symmetry) in the uf sector if we neglect the Yukawa 
couphngs. This approximate symmetry S3 is nearly exact because of the smallness of the 
symmetry breaking Yukawa couplings Y^^. As discussed in the Introduction, the exact global 
Z3 symmetry itself is plausibly assumed to be broken by nonrenormalizable operators such 
that a cosmological domain wall problem does not arise. 

At tree level there is an additional global symmetry in the phase in which the electroweak 
symmetry is unbroken where Hi = and all electromagnetically charged scalars vanish. 
Hence, in the high temperature phase in which the nonsinglet fields are assumed to be 
frozen at their classical potential minimum, we have an enhanced symmetry in the effective 
potential as a function of the singlets only. The enhanced tree level symmetry is Tj-^^Tj^^Tj^, 
where each singlet can be phase rotated independently: 

This symmetry appears because we have tuned the superpotential i^ii^2^3 coupling to vanish. 
Unlike the approximate S3 symmetry, this high-temperature phase classical symmetry has 
significant breaking at 1-loop from perturbative interactions even about the electroweak 
symmetry preserving minima. Nonetheless, this (^3)^ will be useful in understanding the 
SFOPT in which there is a rotation in the singlet sector space during the phase transition.^ 
In addition to the Yukawa and gauge couplings, there are 19 adjustable parameters in this 
model which are \\, n,m\ _^ - = _^_^,m?7 „, Mj, a„ 1^ a kI- Because the neutrino Yukawa 
couplings control the neutrino Dirac mass via the up-type Higgs VEV, these Yukawas are 
small for reasonable values of tan /3, 

and will play a negligible dynamical role. The neutral scalar components belonging to the 
fields {Hi, H2, Li, vf} have finite temperature VEVs denoted by 

<i/°>^ = t;.(T), (uf)^ = v^^AT) J €{1,2,3} (7) 

where the braces represent evaluating the field at the global minimum of the thermal effective 



potential Eq. (16). To maintain the S3 symmetry in the electroweak symmetry breaking 
vacuum at zero temperature and to avoid unnecessary complexities in the SU{2)l charged 
sector, we choose the sneutrino VEVs to be independent of generation, such that 

{v^{0) = real,t;^,(0) = t;^(0) = Tea\,v,c{0) = Vi,c{0) = real}. (8) 

To fix the VEVs in this way, we solve the potential minimization condition for the four param- 
eters {m'\j_^,rn?jj^,m\, m?c}. In addition to the desire to simplify the phase transition history. 



One may also wonder whether omitting the ViV2^^ term is radiatively stable. It turns out that this term 
is generated at 2-loop order, which means that as far as one-loop analysis of this paper is concerned, this 
term can be omitted self-consistently. However, this must be viewed as fine tuning motivated by staying 
consistent with Rcf. [5]. 
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one of our main motivations in choosing f (0) = Vj^c (0) is to preserve the §3 symmetry mani- 
fest in Eqs. ([s]) and Q. Interestingly enough, as we will see, this §3 symmetry spontaneously 
breaks at finite temperature. The left-handed sneutrino VEV is small, as we argue below, 
and the Higgs VEVs satisfy vl{Q) + f ^(0) + 3t;?(0) ^ (0) + t;|(0) = f2(0) = (174 GeVf. 
The rest of the parameter specification will be discussed in Sec. |5| 

To study the electroweak phase transition we need to calculate the thermal effective 
potential as a function of temperature and the field directions which participate in the 
phase transition. There are no charged scalars with VEVs at zero temperature and we assume 
that there are no charge or color breaking minima to appear at finite temperature. In general, 
the left-handed sneutrinos receive VEVs f;>(0) and participate in electroweak symmetry 
breaking, but these VEVs must be much less than the electroweak scale to avoid excessive 
stellar energy loss by v emission |i22j. Hence, we can neglect O {vc) contributions and reduce 
the relevant degrees of freedom to the five dimensional complex field space {ifj*, 5 ^i^}- 
Although a part of the complex phase degrees of freedom in the Higgs sector is a gauge 
degree of freedom, for simplicity we will use the notation {H^, H2, v^} and keep in mind that 
there are nine real degrees of freedom. 

We compute the zero temperature effective potential as a loop expansion over the field 
space {i^i , i?2 > ^i}- The leading order term is the tree level potential given by 

H = ml^ \Hl\' + m% \Hl\' + mj. + ^ {K ' K)' 0) 



0|2 
2 



E 



We exchange the three parameters {rn\j^^m\^^ mlc } for the real VEVs {vi{0),V2{0),Vi^c{0)} 
by solving the three minimization equations 



dVo 



= m\^vi + ^ {yl — t^i + 3 \_—axV2V^,c — KAf2f?c] + 3A^fif2 + 9A^fiw 



VEV 



dm 



= = mjj^V2 - ^ (^1 - vl) f 2 + 3 [-aAfi^i^c - kXv^vIc] + 3\'^V2vl + 9A^f2f?c 
VEV ^ 

= mlcVuc + [-a\ViV2 + (f ?c) - 2KAf if 2t;^c] + 3A^ [vl + vl) Vyc + 2k^VucvIc 

(10) 



VEV 



where "VEV" represents evaluating the fields at the zero temperature vacuum {-ffi , -^2 > ^t} — 
{vi(0), f2(0), f,>c(0)}. Terms proportional to and are negligible and have been omitted. 
Because of the §3 permutation symmetry of our potential, the three equations associated 
with the sneutrino field directions are identical. 

The one-loop radiative correction to the effective potential is given by the Coleman- 
Weinberg potential [23] as a function of the field-dependent mass matrices calculated in 
the Landau gauge ((^ = 0). The mass matrices which we use are included in Appendix [A| along 
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with Hi, the degrees of freedom associated with each matrix that correspond to suppressed 
indices (negative for fermions). Regulating UV divergences in c? = 4 — 2e dimensions, the 
Coleman- Weinberg potential becomes 



^^° = ^E^^TrMf (log 



Ml 



3 
2 



UV 



where Cuv = ^ — 7_e + ln47r and fi is the t'Hooft scale. We impose a mixed renormalization 
scheme in which the counterterms for the parameters {'>TT''jj_^, mj^^, mlc} are chosen such that 
the zero temperature vacuum is unshifted by the radiative corrections. This condition is 
equivalent to requiring tadpole graphs to vanish and imposes 



Srrir 



1 d^V^ 



1 dAV^^ 



V2 d {HIT 
1 dAVl> 



VEV 



VEV 



VEV 



:i2) 



The remaining parameters are determined by the DR scheme and all parameters are specified 
at a renormalization scale of = 100 GeV. We make no assumptions about dominant 



contributions to the one-loop corrections but instead calculate Eq. (11) by summing all 
species that couple to the Higgs sector. 

At finite temperature, the effective potential receives an additional one-loop correction 



2^ 



. b f 



(13) 



where the sums run over bosonic (6) and fermionic (/) mass matrices. The thermal functions 
can be expressed as a a sum of modified Bessel functions of the second kind. 



J2 ^y^^ (^v^) 



(14) 



n=l 

oo 



•Jb {y) = j dxx^ log (l - e 

Jf {y) = dxx' log (1 + e-V-'+y) = - ^^1/^2 (n^) • 

Because these integrals are computationally taxing, we use the Bessel function representation 
and truncate the sum at five terms. This is a very good approximation and introduces less 
than 1% of error. 

At high temperatures, the perturbative expansion fails unless higher order "daisy" graphs 
which diverge quadratically with temperature are resummed. This procedure effectively 
replaces the bosonic field-dependent mass matrix with + where life oc and 
amounts to including a term into the potential given by 



AVr, 



daisy 



T 

12^ 



J2n,Tr {Ml + Ii,f"-{Ml) 



r2\3/2 



(15) 
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The thermal mass corrections life are included in Appendix |Bj Combining all of the radiative 
and finite temperature terms, the one-loop finite temperature effective potential plus daisy 
resummation is given by 

(T) = K) + Al^° + ^V^ (T) + AKiaisy (T) . (16) 

The main difference between the NMSSM and the yUz/SSM relevant for strongly first order 
EWPT can be summarized as follows: 

1. Because of the multidimensionality of the singlet field space {t>;>c(T)}, there can be 
electroweak phase transitions accompanied by rotations within the singlet field space. 
This opens up a new class of phase transitions that are unlike any of the NMSSM 
transitions. For example, the phase transition can take place with the singlet VEV 
hopping from one nonzero value to another: 

{t;,c(T,) = X, ^ 0, {H,) = 0} ^ {v,c{T,) = ^ x„ {H,) ^ 0}. (17) 

The terms in Eqs. ([s]) and ^ that will play a particularly important role for this 
rotational hopping are the soft terms —axHi ■if2^'f + (t'f)^ (which control the cubic 
and lower dimension tree level couplings which in turn control radiative corrections) 
and the superpotential terms —XHi ■ H2i^i + (z^f ) (which control the quartic and 
lower dimension tree level couplings).^ 

2. There is a soft term coupling the singlet to the up-type Higgs, 

AV^oft ~ a^^.Hpl (18) 

which potentially provides a cubic coupling for the Higgs sector. Unfortunately, this 
term does not play an important role in the analysis because Vi, ^ (9(GeV). 

3. The superpotential has a Yukawa coupling of the singlet to the left-handed lepton and 
Higgs, leading to the following additional F-terms: 



AV> = Y,Y, 



K 



(19) 



Given the smallness of the Yukawa couplings ~ (9(10^''), these terms are not 
particularly important for the phase transition when the transition occurs with VEVs 
of order 10^ GeV or less. Note that all of these terms are quartic in nature owing to 
the absence of dimensionful parameters in the superpotential. Also, since the origin of 
R-parity violation is the leptonic Yukawa coupling which is also the source of AVf , 
we see that these do not play a significant role. 

Hence, a generic feature of the /it^SSM SFOPT not reproducible by the NMSSM is the feature 
due to point 1 above. This will be emphasized in the numerical exploration below. We will 
also find one step transitions, which are qualitatively similar to the NMSSM transitions. 



^ As we will see later, the shift of the field origin will generically generate lower dimension couplings from 
higher dimensional couplings. 
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3. QUALITATIVE DESCRIPTION OF THE DESIRED PARAMETRIC REGION 



A novel feature of the /xi/SSM compared to the NMSSM is the transition depicted in 



Eq. (17). In such cases one can shift the origin of the field such that the phase transition of 



interest occurs from the origin. With such shifted coordinates in mind, we define the field 
to be the radial magnitude 



(20) 



for a phase transition controlled by the potential Vr(0) in which the vector of CP-even 
Higgs scalars attains an order parameter change of v. Explicitly, the strength of the phase 
transition is approximately characterized by the SU{2)l breaking \v\/Tc and not (pc/Tc where 
the critical temperature Tc is defined by the condition Vt^{0) ~ Vtc(0c)- 




Figure 1 : A schematic plot of the finite temperature effective potential at the critical temperature 
of Eq. ( 22 ) . The vertical line represents (p = and helps to visualize the effect of i;^ — t- —(j) symmetry 
breaking effect of the cubic term which is responsible for the bump at cj) > for {E > 0, Fna = 0} 



in Eq. (21). 



The finite temperature corrected effective potential of a real scalar field 
temperature will behave approximately as 



near the critical 



Vt(0) 



-M^+ciT^] 



cff. 



* + -Fna(0, 7") + — ^0^+other temperature dependences (21) 



where ci ~ 0{1) constant proportional to coupling constants responsible for the leading mass 
correction and F^^ is the nonanalytic thermal correction contribution that can lead to an 
effective cubic contribution to the potential. Although in the MSSM Fna plays a significant 
role, with a singlet involved such as in the /ii/SSM, need not play a crucial role. Hence, 
we will set Fna = 0. In this section, we neglect "other temperature dependences" in Eq. (21 ). 
Note that Eq. (21) has > even though at T = 0, symmetry is broken when Egg > 



and satisfies a condition specified below. 

Defining where v{Tc) is the degenerate minimum VEV, we find 
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The potential at the critical temperature is depicted in Fig. [TJ The unusual sign of 
2£'g2g/Aeff - M2 stems from our assumption that the symmetry is broken at zero temperature 
due to predominantly the E^.^ term. This situation turns out to be generically beneficial 
for a SFOPT as we explain soon below. The critical temperature Tc is larger if Egg > 



because in that case, the negative contribution from the cubic term in Eq. (21) is enhanced 
for (j) > which means that the quadratic term which is the leading source of positivity (as </> 
approaches (pc from the left) has to be stronger to cancel the stronger negative contribution. 
Since there will be no positive mass squared at the origin during the phase transition in the 
absence of the cubic term, the mass at the origin has to be also larger for increasing Ecs > 0. 
Explicitly, the mass at the origin (which by construction is our starting point of the phase 
transition) is 

dlVT=TM = (23) 
This mass is identical to the mass at = (pc- We can also understand the VEV 

2i?efr 



A, 



cff 



(24) 



which can be heuristically justified by the fact that the broken phase local minimum results 
from a competition between the cubic and the quartic term (which is the dominant source 
of positivity as ^ (j)'^) at the time of critical temperature when the mass term is again 



controlled by Eq. (23). 



Finally, the strength of the SU{2)l breaking in the transition is given by 



(25) 



where 

v{T,) = <f)J{Q) (26) 

and f{^) is a projection cosine onto the Higgs axis. By definition of the SU{2)l breaking 
transition, 0c/(^) < C'(^^(0)). At T = 0, 0(0) = (0)t=o is related to M through 

M = v/3Eeff0(O) - Aeff02(O) (27) 



where Aeff0(O) > 3-E'eff- Equation (25) thus can be rewritten in terms of 0(0) as 



Tc VAeff0(O); ' 



--fm- (28) 



Aoff </'(0) ) V A,ff0(O)^ 

Hence, the strength of the phase transition is controlled mostly by 2 parameters: 



(29) 

Aeff0(O) V Acff 
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Note that since f{^) < 1, this angular projection function can only enhance the phase 
transition in a limited manner. Requiring ^^^^ be real and V"(0(O)) < V{0) result in the 
condition 



< 



E. 



off 



< 



1 



(30) 



Aeff0(O) - 2- 

Therefore, one should keep in mind that although having a nonvanishing E^.^ is good for 
a strong first order phase transition, the enhancement is bounded. Indeed, this bound is 
approximately satisfied by the numerical analysis, and SFOPT points that we find occur 
when 

~ 1 (SI) 

Aeff(/>(0) 2' ^ ^ 



From the derivation of Eq. (28 ), one can see that Eq. (31 ) corresponds to making Tc as small 
as possible during the phase transition. When ^ ^"fft^^ > \ the origin becomes the global 

is defined with respect 



minimum and the symmetry is unbroken. Note also that because 



to the shifted singlet origin in Eq. (20), 0(0) does not correspond to the radial magnitude of 
the scalar field from the original Lagrangian's field origin. 

After this first order phase transition, a second order phase transition might occur when 
V"{Q) = 0. However, with > 0, this does not occur for this ID toy model. Note that 
assumes that there exists a temperature for which V" = which in fact never occurs for this 
toy model. 

Generically, we are interested in a strong first order phase transition characterized by 



V2 



> 1.3 



(32) 



If the asymptotic conditions 



Aeff0(O) 



1/2 and/or Ci/Xes oo are met, the phase 



transition can be arbitrarily strong. However, the following phenomenological constraints 



prevent /constrain an arbitrarily strong transition in Eq. (28): 



1. 



Global minima shifts can prevent the saturation of Ecs/ [Acff0(O)] = 1/2 for a particular 
underlying parametric path. For example, as one approaches i^es/ [Aeff0(O)] = 1/2 
within a particular region of underlying parameter space, ^ the origin of Eq. (20) 



has to be shifted to a new global minimum (where the electroweak symmetry is still 
not broken, i.e. v = 0). When this occurs, {E^s, Xes,Ci} of Eq. (21) undergo a 



discontinuous change as a function of the underlying parameters such as those of 
Eqs. (g and (|4|). 

2. Small Aeff can result in phenomeno logically unacceptably light Higgs (or other scalar 
masses). For example, it is clear from the effective model that 



2A, 



eff^ 



^(0) 



3 E, 



off 



2Aeff0(O) 



(33) 



^ Recall that -Eeff/ [Aoff0(O)] are effective parameters derivable from underlying Lagrangian parametrized for 
example as Eqs. (Isl) and JIl 
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where the term in the parenthesis in Eq. (33) is positive since < -Eefr/ (Aeff(/'(0)) < 



1/2. Note that in any models that embed the MSSM, there is a minimal contribution 
to Aeflf from the D-terms that also makes it difficult to make it arbitrarily small. Note 
also that increasine . lowers the 6 mass as well. 

When _Eeff/(Acfr0(O)) — )■ 1/2, the energy difference AV between the false vacuum and 
true vacuum asymptotically vanishes. Explicitly, we have as Ees Ec = Aefr0(O)/2, 
we find 

AV ^ ^ ^3^f > (34) 



Tc V Acfi / V AEcs 

where AE^.^ = E^s — E^. Since the validity of this estimate requires AV > 0, this 
region of parameter space becomes very sensitive to radiative corrections. 

4. The contributions to Ci that maximize ^/ci/Xcs typically contribute to Xcs as well 
(with different powers). Hence, particularly in the /iz/SSM, we are in a region where 
Aeff is on the larger side and not the small side. 

The features just discussed qualitatively explain the numerical scan of the parameter space 



which identifies a particular parametric region in which Eq. (32) is satisfied at the same 
time with some basic phenomenological constraints which we detail in Sec. [5j There, more 
analytic formulas will be given explaining some of the features of the numerical results. 

Now, let us consider the general path of the electroweak phase transition. At T ^ 
0{TeV), the global minimum will be at 

{v,.{T) = 0,Vi{T) = 0}, (36) 

the scalar field origin.'' As explained previously, the left-handed slepton VEVs are undergoing 
small energy scale transitions which are not particularly relevant to most of the discussion. 
As the temperature is lowered, a nontrivial singlet VEV configuration will realize a global 
minimum, and the system will consequently make a transition. This transition in the sin- 
glets is sometimes accompanied by an electroweak symmetry breaking phase transition and 
sometimes not. If the ffist nontrivial singlet transition is accompanied by strongly first order 
electroweak symmetry breaking, these would be SFOPT from the scalar field origin: 

{v,.iT^) = 0,v,iT+) = 0} ^ Kc(T-) ^ 0,v,iT-) ^ 0}. (37) 



In this case, the origin of the vector whose magnitude is taken in Eq. (20) will be zero. 

In addition, there will generically be singlet transitions from the origin at temperature 
To ffist without an electroweak phase transition, of the form 

K.(T+) = 0,VjiT^) = 0} K(To) = x,(To) ^ 0,VjiTs) = 0}. (38) 



^ Symmetry restoration typically occurs as long as there are no tadpole contribution proportional to the 
temperature. 
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Even if this is a first order phase transition, it will typically complete before the subsequent 
electroweak symmetry breaking, and thus it does not, to leading approximation, participate 
in EWBG. However, it can in principle be relevant for gravity waves (see e.g. [21HSI])- 
Afterwards, there is a subsequent electroweak symmetry breaking phase transition 

{v,.{Tt) = x.(T,) ^ 0, v,{T^) = 0} ^ {vu^iT-) =y.^0, v,{T-) ^ 0} (39) 

whose strength is important for EWBG. In this case, the origin of the vector whose magnitude 
is taken in Eq. (20) will be {vyc{T^) = Xi,Vj{T^) = 0}. When Xi i/i, this transition 



corresponds to a rotation" of the singlet vector. 

Before concluding this section, let us briefly describe how the discrete symmetry discussed 
below Eq. ^ and zero temperature radiative corrections plays a role for some of our strong 



multistep transitions. Once a phase transition of the form Eq. (38) takes place, the set of 
degenerate global minima will form a coset representation of Z3 ® CP ® §3.^ Because of the 
approximate Z3 (g) Z3 (g) Z3 symmetry described in Eq. ([s]) , to tree level accuracy, the coset 
space will be actually bigger: Z3 ® Z3 ® Z3 ® CP ® S3. Some of the Z3 ® Z3 ® Z3 minima 
will be split due to the zero temperature radiative corrections, and the global minimum will 
be at a subset of the Z3 (8> Z3 ® Z3 minima [one of which is what we labeled as x{To) in 



Eq. (38)]. Finally, when the temperature drops enough to make one of the EWSB minima 



degenerate with x(Tc), the transition depicted by Eq. (39) occurs. 

In Sec. [5} we will discuss explicit examples of both one step and multistep phase transi- 
tions. 



4. WEAK SPHALERON AND THE SINGLET 

After the baryon asymmetry has been created at a first order electroweak phase transition, 
it may be washed out by the B-violating sphaleron process [Sj |52] in the broken phase. 
The sphaleron is a nonperturbative field configuration in the Weinberg-Salam theory that 
interpolates between topologically distinct vacua and violates B + L. To avoid washout, 
one must require that sphaleron transitions are suppressed meaning that the rate of these 
processes is less than the Hubble parameter at the time of the phase transition. This imposes 
a lower bound on the sphaleron Euclidean action £'sph(Tc) /T^ > 45 [4j which, in the Standard 
Model, becomes a lower bound on the Higgs VEV in the broken phase \/2v{Tc)/T(. > 1.3 
where v{0) = 174 GeV. The six sneutrino fields of the /xz/SSM which receive VEVs during 
EWSB could in principle modify this bound. As we will see, the modifications are small 
because 1) the left-handed sneutrino VEV is much less than the electroweak scale and 2) the 
singlet sneutrino has a nearly homogenous solution which stays nears the minimum of the 
potential.^ 

To obtain the sphaleron action at finite temperature, we calculate the zero temperature 
sphaleron and apply the scaling law [31] 

i?sph(T)=Esph(0)^ (40) 



^ Recall that §3 is nearly an exact symmetry because of the smallness of the leptonic Yukawa coupling. 
^ This second conclusion was pointed out by Ref. in which the author calculated the weak sphaleron in 

the SM extended by a real singlet. 
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which introduces less than a 10% error. Additionally, we compute the sphaleron solution 
using the tree level scalar potential Vq and neglect radiative corrections. To a very good 
approximation [52] we can also neglect the U(l)^ gauge coupling and compute a purely 
SU(2)^ sphaleron solution. The sphaleron ansatz is static and possesses an SO (3) rotational 
symmetry. The ansatz is given by 



JJCO 



9 

1 / z X + iy 
r \ —X + iy z 



(41) 
(42) 



in terms of the dimensionless radial coordinate ^ = t/tq rescaled by tq = 

(^g^/2^/^^fTvi + vfj ^ {gy/2vy\ All VEVs are evaluated at zero temperature and 

we have introduced = Vy and = Vyc for convenience. We have used the §3 symmetry 
to equate the functions that describe the sneutrino fields of different generations, such that 
the sphaleron solution is given by five functions hi (^) and / (0- With this ansatz the field 
equations become 



vf dhj 



eh'i+2^h',-2{i-ffh, 

3 leK + '^iK] 



1,2 

3 

4 



1 



er - 2/(1 - /)(i - 2/) = --e{i - /) [vihi + vihi + ^vi^^ g\i m 

where the potential is normalized to vanish in the EWSB vacuum Vq ihi = 1) = 0. Note that 
the term (1 — /)^/ij is absent for z = 4 because the singlet sneutrinos do not couple to the 
gauge bosons. The sphaleron action is obtained by integrating the sphaleron solution 



£^sph (0) 



47rt;v^ 



d[_ 
di 



+ 



2v' 



dhx 
d( 



+ V2 



dh2 



+ ^ {vfhl + vlhl + ^vll^^ (1 - ff di 







4 


2 


^ a 




V3- 


-If 





d^ 



+ 3 U4 



dh^ 



(44) 



We can observe immediately that contributions from the left-handed sneutrinos will be neg- 
ligible because the function always appears with a prefactor of = Vc v . 



We can study the sphaleron solution by considering the asymptotic limits of Eq. (43). 
In the large ^ limit, all five field profiles must asymptote to unity in order for Egpi^ to be 
finite. In the small ^ limit, we find that the gauge boson and the three weakly charged scalar 
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functions asymptote to zero, as in the Weinberg-Salam model, but that the singlet function 
approaches a value which can in general be nonzero: 

h^iO^M ^e{l,2,3}, h,^c + p£. (45) 

The singlet function behaves differently in this limit because the gauge coupling term (1 — 
/)^/i4 in the field equation is absent. The boundary condition on the singlets makes the 
solution for /i4 (^) qualitatively different than for the Higgs fields. In particular, the solution 
/i4 (^) which minimizes i^sph will tend to be homogenous with h4 ^ 1 for all ^. The solution 
is homogenous because E^p^ 3 {dh4/d^ f is positive semidefinite. Hence, it can be minimized 
by a constant /i4, and the solution remains near = 1 because this is where Vq is minimized. 
As a result, the singlet fields contribute negligibly to the sphaleron action. 

The sphaleron solution and energy density for a fiducial parameter set are plotted in Fig. 



2l To obtain the field profiles we solve Eq. (43) in the large and small ^ limits analytically, 



then match the solutions at five radii which are chosen to minimize -Esph- As discussed 
above, the singlet solution hovers around = 1 where the potential has a minimum. To 



display how each term in Eq. (44) contributes to the sphaleron action, we have also plotted 
the integrand for the gauge kinetic, scalar kinetic, and scalar potential contributions sepa- 
rately. We observe that the sphaleron action is dominated by the kinetic terms. Since the 
parametric dependence only appears explicitly in the scalar potential, which is negligible, 
we expect that the sphaleron action is largely independent of our parameter choice. For this 
parameter set we find £^sph(0) ~ 1.83^ ~ 8.7 TeV which translates into a bound on the 

Higgs VEV at the critical temperature that is > 1.3. As such, the Higgs VEV must 

satisfy the same constraint in the /iz/SSM as in the SM to avoid washout. 



£,pl,(0) / (4;rvV2"/g) 
1.0 




Figure 2: On the left, the //i^SSM sphaleron solution versus the dimensionless radial coordinate 
with hi and /i2 (solid line), (dashed line), /14 (dotted line), and / (dashed-dotted line). The 
solution /i£;c for the singlet sneutrinos does not satisfy the same boundary condition at ^ — )• as 
the SU(2)^ charged scalars. Hence, the solution of minimum energy is the one in which /ij>c 1 for 
all ^. On the right, the sphaleron energy density, Eq. ( |44[ ), with gauge kinetic terms (dashed line), 
scalar kinetic terms (dotted line) , scalar potential terms (dashed-dotted line) , and the total energy 
density (solid line). This plot illustrates that the sphaleron action is dominated by the kinetic 
terms and that the contribution from the scalar potential is negligible. 
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5. PARAMETER SCAN AND PHENOMENOLOGICAL BOUNDS 

We have investigated the /xz/SSM phase transition by performing a two-dimensional pa- 
rameter space scan. For the two free parameters we use rrich = SXvj^c, which coincides with 
the charged Higgsino mass in the limit M2 ^ miy,mch, and the dimensionless variable 

These parameters are scanned uniformly over the ranges rrich '■ [75 GeV, 175 GeV] 
and a : [0, 25] by varying Vj^c and ax- The SUSY-breaking parameters are cho- 
sen to match the conventions of [Sj: |m^, m|c, m|^, mgc| are fixed at a fiducial 

SUSY-breaking scale which is taken to be 1 TeV, gaugino masses are set to 6M1 = 
3M2 = M3 = 3 TeV, and A-terms are scaled by the associated Yukawa couplings as 
{Au^ = au/Yu^,Ad^ = ad/Yd^,Ae^ = ae/Y^^.A^ = -a^/Y^} and fixed at 1 TeV. Note that we 
have assumed for simplicity common left-handed sneutrino VEVs f ,7. = Vy and a diago- 
nal Yukawa matrix {Yy)^- = Y^Sij. The remaining soft masses {m|^^, m^^, m|, m?c} are 
exchanged for the VEVs {vi (0) ,V2 (0) ,f£; (0) ,f£;c (0)} by solving the minimization equa- 
tions at zero temperature. The remaining Higgs sector parameters are chosen to be 
tan P = 2.6, k = -0.64, A = 0.18, t;^ = 1.4 x 10"^ GeV, and a« = -236 GeV. Given that some 
of our sparticle masses are far larger than ~ C(100)GeV, we could have integrated out 
these fields giving rise to a more illuminating effective field theory parametrization within 
the DR scheme.^ However, to stay similar to the parametrization used in [UlS], and to give 
a relatively unrestricted range for possible Tc, we have kept these relatively heavy fields as 
dynamical. 

At each point in the parameter space, we calculate the /xz/SSM spectrum. In order 
to get a handle on phenomenological constraints, we impose the MSSM search bounds 
for the SUSY particles and require the Higgs masses to be > 90 GeV [HH] (later we will 
show a sample parametric point Higgs spectrum with the lightest Higgs mass of about 110 
GeV). Model dependent bounds are of interest, but typically they are weaker, as far as 
the neutral Higgs is concerned, because of singlet mixing effects. A more complete model 
dependent phenomenological consistency check including the study of charged Higgs me- 
diated 6 — > 57 rates is beyond the scope of this paper. We calculate the spectrum of 
the charged Higgses (0f ), charginos {xt)^ and neutralinos (Xi) at tree level and require 

|m^± > 79.3 GeV, m-± > 94 GeV, m-j^o > 46 Gev|. The SM-like neutrinos mix with the 

neutralinos and heavy neutrinos in a seesaw matrix. We are able to reproduce the correct 
neutrino mass scale but neglected the question of precise neutrino mass pattern*^ since any 
desired neutrino mass pattern will not be difficult to achieve by adjusting the small Yukawa 
couplings. Since we have already noted that the smallness of the leptonic Yukawa couplings 
make their role in the current SFOPT analysis insignificant, this does not present a significant 
loss of generality. 



^ Recall that in DR scheme, decoupling is accomplished "by hand" through computing threshold corrections 

after integrating out fields. 
® The issue of neutrino masses in the ^i^SSM was studied more extensively in Ref. [561 157) . 
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Because the squark, charged slepton, and left-handed sneutrino masses are supported 
by their TeV-scale SUSY-breaking mass parameters, they are insensitive to parameters in 
the Higgs sector and are not affected by the phenomenological lower bounds. We compute 
the mass spectrum of Higges and singlet sneutrinos at one-loop order using the effective 
potential. Since we choose the VEVs for these fields to be real and there is no explicit CP- 
violation, the CP-even (0j) and CP-odd (a,) components do not mix. The mass matrices are 
given by the curvature of the one-loop effective potential evaluated at the zero temperature 



vacuum 



(Mi). 



(92T/0 



d{TZipi)d{TZ(pj 



[Ml) 



VEV 



d{Xipi)d{Xipj 



(47) 



VEV 



where y^^ G {HI,HI,p'^]. We can separately impose the mass bounds m^^ > 92.8 GeV and 



ma, > 93.4 GeV. 



At each point in parameter space that satisfies the phenomenological mass bounds, we 
require the electroweak-breaking vacuum with v{0) = 174 GeV to be the global minimum of 
the one-loop effective potential. This condition imposes particularly strong constraints on 
the parameter space. To understand these constraints and the nature of our multistep phase 
transitions, we must discuss the structure of the {i^?, H2, z>f } field space and, in particular, 
determine the locations of the vacua that could potentially have lower energy than the EWSB 
vacuum. Recall that in the subspace with = = there is a Z3 ® ® Z3 symmetry 
at tree level given by Eq. ([s]). To locate the extrema in this field space we solve the three 
cubic equations 



duf 



ze {1,2,3} 



(48) 



HO=hO=o 



for z/f. The solutions of Eq. (48) which turn out to be minima (for our choice of the sign of 
Ok) are given by 



I'i = pie 
Pi = or 



G {0,1,2} 

1 

4^ 



(49) 



We will focus on the solutions with pi = p2 = Ps ~ v^c because these minima are in general 
deeper than those with pi = 0. Then, there are 3^ = 27 local minima in the = H2 = 
field space, that we will refer to by Xmn2n3 where the subscript indicates the phases of the 
three singlets. The (^3)^ ® §3 ® CP symmetry ensures the degeneracy of the 27 minima. 
Moving away from = H2 = as illustrated in Fig. [3, the (^3)^ symmetry is broken to 
Z3 by terms in Vq proportional to a\ and A. We will use ynin2n3 to denote the point in field 
space near to x„j„2n3 but where H^/vi = H2/V2 = 1. For example, in this notation, yooo 
corresponds to the EWSB vacuum in which the three singlets have real VEVs. 



^ Since the effective potential is defined as a sum over IPI diagrams with zero external momentum, this 
definition of mass differs from the pole in the propagator by the difference of the scalar self-energy evaluated 
at p ^ rripaic and p — 0. 
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Figure 3: The tree level potential plotted over a slice of the field space with = H2 = 
on the left and H^/vi = H2/V2 = 1 on the right. The labeled points are defined in the text, 
and a stationary point of the potential can be found at or near each of the labeled points. The 
potential grows farther from the central region. In the EW-preserving subspace, the three minima 
are degenerate, but the Higgs VEV selects out yooo as the global minimum. 

At one-loop order, radiative corrections break the approximate (^3)^ symmetry [described 
above Eq. ([s])] and split the degeneracy of the Xn^rij^g minima as represented by Fig. |4j After 
including radiative corrections, the preserved symmetry group is Z3 S3 ® CP. (Here, as 
an approximation, we are ignoring the fact that the subgroup §3 ® CP is explicitly weakly 
broken in the Lagrangian already while Z3 must be broken by nonrenormalizable operators 
to evade cosmo logical inconsistencies caused by domain walls.) The 27-fold degeneracy is 
split into three classes: a 3- fold degeneracy of the points Xm, a 6-fold degeneracy of the 
points Xijk for i ^ j ^ k, and a 18-fold degeneracy of the points Xuj plus permutations for 
i ^ j- In order to discuss the phase transition we will choose one representative from each 
class: afooo) 2:0125 and Xqoi- In this notation, if we say a transition occurs from the origin to 
X012 we mean that just below the critical temperature the vacuum is localized nearby to one 
of the six field points in the class that contains xoi2- 

The radiative corrections will generally split the degeneracy in such a way that some of 
the EW-preserving vacua will be depressed relative to the EWSB vacuum and may cause 
the latter to become metastable. This is both good and bad for the parameter space scan 
and phase transition analysis. It is bad because many points will be excluded because the 
EWSB vacuum is only metastable. On the other hand, it is good because with appropriate 
tuning, we can obtain an EW-preserving vacuum that is nearly degenerate with, but slightly 
higher than, the EWSB vacuum. Along the trajectory connecting these vacua, we can 
make E^^/ (Aoff0(O)) arbitrarily close to one half and obtain SFOPT. In Appendix [C| we 
include analytic bounds which must be satisfied to prevent the EWSB vacuum from becoming 
metastable. 

At each point in parameter space which satisfies the mass and vacuum bounds described 
above, we calculate the critical temperature, Tc, and Higgs VEV, f(Tc), at the electroweak 
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Figure 4: Same as Fig. [3] but the contours represent values of the one-loop effective potential 
at zero temperature. The degeneracy is broken even in the EW-preserving subspace and induces 

n° (^012) < ^1° (xoOl) < 14° (^OOO) . 

phase transition. The phase transition is calculated using the following procedure: increase 
the temperature from zero in increments, at each temperature minimize the thermal effective 
potential to find the EWSB vacuum VEWssiT) = [vi{T),V2{T),Vi^^(T)^ , also find whatever 
metastable vacua VMS,i are near to Xn^n2nz and ymn-iuzi as the temperature increases the 
location and depth of these stationary points will change, converge on the critical temper- 
ature Tc at which the EWSB vacuum becomes degenerate with one of the EW-preserving 

minima Vj^ {vewsb{Tc)) = Kff {vMS,iATc)), compute v (TJ = ^J\vl{T^)f + \v2{Tc)f as the 
Higgs VEV in the broken phase. Using this procedure, we obtain Tc and v (Tc) for the lowest 
temperature phase transition. Generally in this region of parameter space, multiple phase 
transition steps are required to bring the field configuration from the high-temperature sym- 
metric phase to the zero temperature broken phase. We must investigate separately earlier 
steps. 

The results of the 2000 point parameter space scan are summarized in Fig. [5]where regions 
Ilia and Illb are the only likely viable regions for SFOPT EWBG. We will describe the 
different regions here and give an analytic derivation of the boundaries and their parametric 
dependence in Appendix [C| The points in region I are excluded because the EWSB vacuum, 
where v{0) = 174 GeV, contains a tachyonic direction. The points in region II are excluded 
because the EWSB vacuum is metastable. For regions Ila, lib, and lie, the actual vacuum 
can be found at the following points: the origin of field space in region Ila, nearby to Xoi2 in 
region lib, and nearby to ?/oi2 in region lie. That is, in regions Ila and lib, the electroweak 
phase transition does not occur. Region lie does not work for EWBG as well as we will see 
below. In region III there are no tachyons, no false minima, and all phenomenological mass 
bounds are satisfied, but as we will see only Ilia and Illb are likely to give acceptable phase 
transitions for EWBG. 

The phase transition at each point can be classified into one of four types based on the 
path that the vacuum follows through the {H^, H2, i^f} field space. In the largest region 



20 




Figure 5: A slice of the ^z^SSM parameter space. Region I suffers from tachyons in the EWSB 
vacuum. In region II the EWSB vacuum is metastable and we exclude these points. In region III 
we calculate the electroweak phase transition and find that the path through field space can be 
classified into one of four types, shown on the right. 



Ilia, the PT makes two steps: from the origin to a |EW, ^3, ^3, Q't'} phase and then to 
the {eM^,^3,§3,CP} phase. In region Illb the EWPT occurs in one step directly from the 
origin to the EW-broken phase. In region IIIc, the EW symmetry is broken by a second 
order phase transition in which only H2 gets a VEV; then, a first order phase transition 
occurs giving the singlets VEVs. Finally in region Illd the phase transitions occur in three 
or four steps and there are multiple EWSB phases, whose details for a representative point 
are discussed below. However, as we will see, region Illd is unlikely to give an acceptable of 
EWBG scenario. 

To understand how the /iz/SSM phase transition differs from the NMSSM scenario, we 
have taken one representative parameter point from each sector of region III and followed the 
full phase transition from the origin xq to the zero temperature EWSB vacuum ^000 • In the 
tables, the minima above and below an arrow are degenerate at the temperature indicated. 
A 0^ indicates that a second order phase transition occurs along the specified field direction. 



Ilia. Two Step Transition via EW-preserving Phase: xq > ^012 > Vooo 
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Representative point: {mch, cr} = {108.8 GeV, 9.12}. 



T 


{Hf,Hli)f} (GeV) 


V2'-P 


Symmetries 


> My, 

75.1 GoV 
> 


{0,0,0,0,0} 
{0, 0, 191.7, 191.7e^2V3^ 191.7e*^''/3| 






{EW,Z3,S3,CP} 
{EW,%,^3,Qf} 


54.6 GeV 

y 


{O, 0, 192.6, 192. ee*^'^/^, 192.6e*'*''/3} 
{61.7, 160.7, 201.0, 201.0, 201.0} 



4.46 


{Ew,%,^3,gf} 
{eM^,%,S3,cp} 


GeV 


{62.5, 162.5, 201.5, 201.5, 201.5} 


N/A 


{eM^,%,S3,cp} 



Table I: Phase transition path for the representative point in region Ilia. 



At T = 75.1 GeV, a first order phase transition gives the singlets VEVs and breaks 
Z3,S3, and CP. The EW symmetry is broken by a strongly first order phase transi- 
tion at 54.6 GeV which also restores S3 and CP. Baryon number may be generated 
at the strongly first order EW-breaking PT because sphalerons are suppressed by 
y/2v{T^)/T^ = 4.46 inside the bubble. 

-, IPT 

Illb. One Step: Xq > Vooo 

Representative point: {rrich, o"} = {102.5 GeV, 7.22}. 



T 




V2^-^ 


Symmetries 


> 

69.6 GeV 
> 


{0,0,0,0,0} 
{60.3, 157.7, 188.3, 188.3, 188.3} 



3.43 


{EW,Z3,S3,CP} 

{eM^,%,§3,cp} 


GeV 


{62.5, 162.5, 189.8, 189.8, 189.8} 


N/A 


{eM^,%,§3,cp} 



Table II: Phase transition path for the representative point in region Illb. 



At T = 69.6 GeV, the Higgs and singlets obtain VEVs simultaneously breaking the 
EW symmetry and Z3. This one step phase transition resembles the ones seen in 
certain parametric regions of the NMSSM and other Higgs-singlet extensions. A baryon 
number may be generated since y/2v{Tc) /Tc = 3.43 in the broken phase will suppress 
washout. For the parameters in region Illb, we have plotted in Fig. [6] the order 
parameter and critical temperature as functions of -Eefi/ {^es4'{0)) which we calculate 
using the tree level potential along the trajectory that joins xq and yooo- The order 
parameter grows and the critical temperature decreases as -E'es/ (Aeff0(O)) approaches 
1/2 from below. The data points do not extend all the way to 1/2 because the radiative 
corrections lift the potential in such a way that parameter sets with Ees/ (AcfT0(O)) ~ 
1/2 at tree level have a metastable EWSB vacuum at one-loop. 

2PT IPT -, 

IIIc. Two Step via EWSB Phase: xq > Vh-j ^ ^000 

Representative point: {mch, cr} = {95.1 GeV, 5.37}. 



At a high temperature T = 116.4 GeV the EW symmetry is broken by a second order 
phase transition along the up-type Higgs direction. As the temperature decreases, the 
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T 


, J. J. 2 ) "j J 


V2^ 


Symmetries 


> Mw 

116.4 GeV 
> 


(0,0,0,0,01 
{0,0+, 0,0,0} 




0+ 


{EW,Zs,Ss,CP| 
{eM^,Z3,S3,CP} 


57.6 GeV 

> 


{0, 101,0, 0, 0} 
{61.6, 160.6, 175.4, 175.4, 175.4} 


2.5 
4.25 


{E/V,%,§3,CP} 

{eM^,%,S3,cp} 


GeV 


{62.5, 162.5, 176.2, 176.2, 176.2} 


N/A 


{eM^,^3,S3,cp} 



Table III: Phase transition path for the representative point in region IIIc. 



global minimum of the effective potential moves along the H2 axis until it becomes 
degenerate with a minimum localized near to |/ooo- A first order phase transition occurs 
with \/2v{Tc)/Tc = 4.25 inside the bubble and \/2v{Tc)jT(. = 2.5 outside the bubble. 
In this scenario, there is no baryon number generation. Because the first transition 
is of the second order, there is no coexistence of phases. The second transition is 
first order, but the sphaleron transition rate is suppressed both inside and outside the 
bubble such that B + L is preserved on both sides. 

Illd. Multi-Step: xq > ^012 > {Vooi or ^002) > Vooo 

Representative point: {mch, c} = {121.6 GeV, 4.40}. 



T 




V2^ 


Symmetries 


> 

222.3 GeV 
> 


{0,0,0,0,0} 
{0,0, 197, 197e^2V3^ 197e^'^''/''^} 






{EW,Z3,S3,CP} 


128 GeV 
> 


{0, 0, 229, 229e^2'^/^ 229e''^''/^} 
{0+, 0+, 229, 229e*2V3^ 229e''^''/^} 




0+ 


{EW,^3,^3,Qi'} 


105 GeV 
> 


{0+, 82, 232, 232e^2V3^ 232e*'^''/-'^} 
|29eii-9'r, 117e*^-9'', 229, 229, 2326^2^3} 


1.1 

1.6 




89.0 GeV 
> 


{32.0, 128.0, 230.3, 230.3, 232.76^2'^/^} 
{58.4, 152.3, 224.4, 224.4, 224.4} 


2.10 
2.59 


{Ep/,^3,f^3,0P} 
{eM^,^3,S3,CP} 


GeV 


{62.5, 162.5, 225.1, 225.1, 225.1} 


N/A 


{E/V,^3,S3,CP} 



Table IV: Phase transition path for the representative point in region Illd. 



At this parametric point, the phase transition occurs in four steps with the EW sym- 
metry broken in the second step by a second order phase transition. As the tempera- 
ture drops from 128 GeV to 105 GeV, the sphaleron becomes increasingly suppressed. 
When the IPT occurs at 105 GeV, the sphaleron is inactive, such that there will be no 
B-number generation. Not every phase transition in region Illd follows this particular 
PT path, but the PTs are generally multistep with at least one EWSB intermediate 
phase and transitional CP violation. 

To give an impression of the particle masses in this region of parameter space, we 
include here the spectrum for the representative point in region Illb where {rrich, c} = 
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Figure 6: The order parameter (squares) and critical temperature (circles) plotted against 
-E'eff/AcfT'/* (0) for the points in parametric region Illb. We calculate £'cfr/-^cff0 (0) using the tree 
level potential along the trajectory that connects the origin xq and zero temperature vacuum yboo- 



{102.5 GeV, 7.22}. The slepton, squark, gaugino, and left-handed sneutrino masses are all 
O (TeV) because we have fixed the soft masses in these sectors at a fiducial SUSY-breaking 
scale. In the case of the z/, we solve for the soft mass using the minimization equations to 
find m| ^ —auVi,cV2/vy, yet it is still typically TeV scale because Ai, = ay/Y„ = —1 TeV. 
The remaining fermion and charged scalar spectrum is given by 

mH± = 342 GeV 

= 96 GeV 
m^c = 260, 243, 243 GeV 
^Xi2 ~ 88.7,94.5 GeV, mostly Higgsino 

= 41,41,7.8 meV (50) 

which are calculated at tree level. The LSP is a Higgsino with mass 88.7 GeV. The degen- 
eracies present in the neutrino sector result from the S3 symmetry of our Lagrangian. By 
allowing the left-handed sneutrinos to have different VEVs or choosing different values for the 
Yukawas, we could obtain a correct neutrino hierarchy. We include these masses here to 
demonstrate that the seesaw matrix produces the correct mass scale for the light neutrinos. 
The neutral scalar masses are calculated at one-loop using the effective potential. Because 
there is significant mixing, we have included their mass eigenvalues and field composition in 
Table |5| Once again the degeneracies are a result of our S3 symmetry in the singlet sector. 
The lightest Higgs is mostly up- type with a mass of 110 GeV at this parametric point and 
only varies by 10 GeV over all of region III. 



6. DOMAIN WALLS 

It is well known [T71 [581462] that domain wall formation can be cosmologically problematic 
when spontaneous breaking of discrete symmetry occurs. In our scenario, we have only one 
"exact" discrete symmetry Z3 at the level of explicit parametrization of the Lagrangian. 
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Mass (GeV) 








Kz>2" 




one; 


0.83 


U.lo 


u.oo 


0.00 


n r\r\ 

o.ou 


140 


0.01 


0.21 


0.26 


0.26 


0.26 


128 


0.00 


0.00 


0.58 


0.01 


0.41 


128 


0.00 


0.00 


0.09 


0.65 


0.07 


110 


0.16 


0.63 


0.07 


0.07 


0.07 



Mass (GeV) 




9^2° 






9z>| 


439 


0.48 


0.07 


0.15 


0.15 


0.15 


369 


0.00 


0.00 


0.02 


0.39 


0.59 


369 


0.00 


0.00 


0.65 


0.27 


0.08 


314 


0.39 


0.06 


0.18 


0.18 


0.18 



Table V: CP-even and CP-odd Higgs masses and mixings for a sample parameter point. The field 
composition is described by the squared eigenvector associated with each eigenvalue. 



Because of the undesirable cosmological consequences of domain walls, we have implicitly 
assumed that this symmetry is broken by nonrenormalizable operators which are cutoff by 
a scale larger than many TeV (such as not to disrupt the effective potential analysis). In 
addition, we have approximate discrete symmetries such as §3 ® CP which a priori can cause 
problems if the symmetry breaking operators are overly suppressed. However, the set of 
electroweak symmetry breaking vacua of interest in this paper does not break §3 ® CP (i.e. 
our symmetry breaking pattern can naturally select a S3 (8> CP singlet VEV to be the lowest 
energy vacuum as partly demonstrated in Appendix D). Hence, we will neglect any transient 
behavior and focus on Z3 domain walls even though the analysis is not very specific to the 
discrete group. Although a full analysis of domain wall histories is beyond the scope of this 
paper, here we briefly estimate the effects of the suppressed symmetry breaking operators 
that will alleviate the cosmological problems associated with domain walls that may form 
when the discrete symmetries considered in this paper are spontaneously broken. We will 
follow closely the work of Ref . |17j . 

In Ref. [17] , it is estimated that during the approximate discrete symmetry breaking phase 
transition, domain walls separating approximately degenerate minima are formed. Then a 
simplified model of domain wall evolution is considered which approximately accounts for 
the surface tension of the bubble, the friction coming from bubble wall interaction with the 
plasma, and the pressure coming from energy density difference between the approximately 
degenerate minima. This last ingredient (pressure from energy density difference) is what 
will be coming from the inclusion of suppressed symmetry breaking operators, and we will 
refer to this simply as "pressure difference." If the pressure difference dominates, one of 
the approximately degenerate phases will eat away at the higher energy phase regions and 
eventually dominate in a time scale controlled by the strength of the symmetry breaking 
operator. 

Estimating the friction to be negligible, an approximate sufficient condition for curing the 
possible domain wall problem from a cosmological perspective is to have the pressure differ- 



Because we have + < 0.5 in the parametric reginie of interest, the couplings should remain pertur- 
bative up to close to the GUT scale [2T] (see e.g. [S] for explicit plots which suggest that our parametric 
choice is close to the border of perturbativity up to the GUT scale). Thus we are not severely restricted 
in the cutoff scale of our effective field theory. 

For example, if one wanted to analyze transient domain walls associated with S3 breaking, one can easily 
work out from our Lagrangian the leading effective scalar operator breaking S3 and use the result at the 
end of this section. 



25 



ence dominate before the equilibrium initial condition period of big bang nucleosynthesis: 
i.e. before the photon temperature reaches about 10 MeV. Explicitly, assuming order unity 
Lorentz factor 7 for the bubble wall speed, one must require 

' > if) 1^1' 

where e is the energy density difference coming from suppressed symmetry breaking opera- 
tors, a is the energy per unit area of the bubble wall, and R{t) is the time dependent radius 
of a typical bubble. For a dimension 4 + m nonderivative operator consisting of scalars only, 
e can be estimated as 

' - (52) 

where A is the cutoff scale and we have assumed all scalar VEVs to be of common order 
V = 174 GeV (which is appropriate for our scenario). To be able to treat m = 0, we will 
set A = 100 TeV and find a bound on the value of Cu for different values of u. Assuming 
R t ^ 1/H (where H is the Hubble expansion rate) and cr ~ v^, we find 

Cu > 10-2^500" (53) 

for u > O}^ Hence, as long as the cutoff is not required to be very large (in contrast with 
the assumption of Ref. [IT]) or the accidental symmetry arising from the UV completion 
quantum numbers do not make u too large, this bound is very easy to satisfy for the Z3 
domain wall problem. Of course, if the cutoff is taken to be high and/or a UV completion is 
desired without fine tuning, model building challenges along the lines of Refs. [SH |M] exist. 



7. SUMMARY 

We have uncovered a /ii^SSM parametric region giving rise to a first order phase tran- 
sition sufficiently strong to be useful for the electroweak baryogenesis scenarios involving 
electroweak symmetry breaking bubbles as the source of out of equilibrium and SU{2)l FF 
operators as a source of baryon number violation. The parametric region corresponds to 



tuning the soft terms in the Lagrangian axHi-H2i'f and — m?c|z>fp to achieve Eq. (31). The 
numerical values of the uncovered parametric region is in the paragraph containing Eq. (46) 
and regions Ilia and Illb depicted in Fig. |5| As expected, the Yukawa coupling of the singlets 
to the leptonic sector does not play a role in determining the strength of the phase transitions 
because of the weakness of the coupling tied to the smallness of the neutrino masses. 

The region Ilia transitions are two-step transitions in which the electroweak symmetry 
breaking is the second transition that starts from a phase in which the singlet scalars of 
the /ii/SSM have a nonzero vacuum expectation value (e.g. starts from a vacuum which 
spontaneously breaks the approximate 5*3 symmetry in the singlet sector). These transitions 
contain a rotation in the singlet field space and do not have an analog in the NMSSM 
transitions because of the different dimensionality in the singlet complex vector space. The 
region Illb transitions are the ones in which electroweak symmetry breaking transition starts 



This result can easily be checked to be consistent with Ref. [17] for u = 1. 
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from the origin of the scalar field space. All these transitions have useful descriptions in terms 
of the representations of the approximate discrete symmetries in the system. 

Our phenomenological bounds were rather minimal and placed using Ref. [S3], but in 
many parametric regions, the observables are sufficiently far away from the bounds that 
the plausibility of the phenomenological self-consistency is strong. Follow-up possibilities 
include a more complete collider related phenomenological investigation in this parametric 
regime, studies of domain wall histories due to weak global symmetry breaking operators, 
and a complete computation of CP asymmetry creation and transport leading to baryon 
asymmetry. 

Given that the /iz/SSM had to give up the popular thermal leptogenesis scenario due to 
its low scale implementation of the type I seesaw, this work is of interest as it shows that 
electroweak baryogenesis may be a promising avenue to create baryon asymmetry in this 
class of models. Given that the /xz/SSM is one of the few supersymmetric models in which 
all dynamical degrees of freedom responsible for the neutrino mass may be accessible at 
TeV-scale colliders, it is encouraging that the model has a good chance at being consistent 
with the observed baryon asymmetry in the universe. 
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Here we include the tree level, field-dependent mass matrices which are required to 
compute the one-loop radiative corrections. We fix the charged scalars at their vanishing 
VEVs, let the left-handed sneutrino VEV z/j = Vy be real, and treat the matrices as functions 
of the complex fields {H^, H2, uf}. 



Neutralinos. Using the basis ix^) = W^s, if^, ifg, J^i, ^^2' ^^35 '^i' '^s, ^'sj, the mass 
term appears as £ 9 — | (x^)'^ M^o -|- h.c. with n^o = —2 and 
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Appendix A: Field-Dependent Mass Matrices 





(Al) 




T 



(A2) 
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and Al is a symmetric, sparse matrix with nonzero elements 



Mm? 



Mb 



MwsWa 



M 



M 



W3H? 



Ml 

92 



Mh«hO 

MffO^c 
MfjO^c 



-X {VI + v^^ + vl) 
-XHl 



(A3) 
(A4) 

(A5) 

(A6) 
(A7) 

(A8) 

(A9) 
(AlO) 
(All) 
(A12) 



Charginos. Using the basis {^^)^ — |— iA+, 
|— iA~, , e£, T^l the mass term appears as £ 9 —^{^~)^M^±{^'^) + transpose 



Ht.^R.^^R.r^] and (*-)^ 



with = 



-2 and 



/ M2 


92 (HIT 








\ 


92 {Hir A 


[Dl + D^^ + DD 




-YiiVi, 


-Y^Vy 


92Vi:, 












g2Vi> 







Y 





\ 92Vi> 










Y^El ) 



(A13) 



Gauge Bosons. The propagators and field-dependent masses in the gauge sector have 
gauge dependence. We work in the Landau gauge (^ = 0), in which the scalar component and 
ghost propagators have no field dependence. The charged gauge bosons have field-dependent 
mass 



9l 



+ |i/°|%3v. 



(A14) 



and in the basis {W3, B) the neutral gauge bosons have the mass matrix 



9l 



9192 
2 




(A15) 



In order to count the degrees of freedom in the gauge sector, we must distinguish lon- 
gitudinal and transverse components of the gauge boson fields, 2n 



n 



W4. 



4 and 

2nw-iBL = ^WsBt — 2. We do this because only the longitudinal components receive thermal 
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mass corrections in the computation of the daisy correction Eq. (15), [65]. Note that vf, is 
numerically negligible in these equations over all parameter regions of interest. 

Squarks. Because we assume there is no intergenerational mixing in the squark sector, 
the squark mass matrix block diagonalizes. The i^^ generation up- and down-type squarks 
have mass terms C 3 —^qjM^.qi in the basis qi = {qLi, Qr-} with rig = 12 and 

{Ml)u = ^i + l{^-§) {K - K + 3-^) + K \H',\' (A16) 

(MDu = i^Dli = (^2)* - Y.Y^v, {vl + z>2^ + vl) - (if?) {vl + v^, + v^) (A17) 

(M^J,, = <^ + 'i (l^?r - \< + H) + \Hl\' (A18) 

H),, = MJ - 1 (M + 1^ (|^o|2 _ |^o|2 + 3.^) + Yl \Hl\' (A19) 

= H)l = " ^'^^ + + ^^^^ 

(^i),,=-!^-f (i^?r-i^°r+3-o+^^-i^?r • (A2i) 

Charged Scalars. The charged Higgs mixes with the charged sleptons. Using the basis 
= [H^ , H2 ,ej^,e~l^, fii, fi~^,fl ,f^^ , the mass term is £ 9 —S^M'^^S" with nH± = 2 
and the elements of the Hermitian mass matrix are 



0|2 
2 



i=l 

+ X^\vl + v^2 + vl\^ (A22) 

= «a E + 3^.1;A [Hiy + (^^,2 - SA^) [H'.Hiy + K\Y, {Vff* (A23) 

= < + ^ isl + gl) |i^,t - \ ial - 91) {H + + E 

i=l 

+ \9l + 9^ + u'sf (A24) 

{Ml^)H,i,^ = -aev, - Y,Y,^ [Hpty (A26) 

{MH^)H,h^ = (^^2 - n') ^.i^2° + Ay.ii?ii° - y.^ iKf - du a^tr im 
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k=l 



,22 



m 



5i 



fc=i 



0|2 



YeYe^i 



[vir 



(A28) 

2 

(A29) 
(A30) 

(A31) 



Neutral Scalars. The neutral Higgses mix with the left- and right-handed sneutrinos 
in a 16 X 16 matrix M\^. At the EWSB vacuum which respects CP, this matrix block 
diagonalizes into CP-even and CP-odd sectors. In order to study the phase transition in 
which there are transitional CP-violating phases, we must retain the off-diagonal blocks. In 
the basis 0^ = {Sfti^o; ^-f^o} where i^o = H2, i^i, 1^2; ^2, ^'3} the mass term is given 

by £ 9 — |0"^M^p0 with uhq = 1- One can obtain the mass matrix 



by differentiating the full scalar potential 



(A32) 



V, = V, + ml'£\h? + Y^\\Hl\'Y. 



V^aly.yt + Y^K {H^,i),y {utf + h.c] 



+ 



gl + gl 



0|2 



E 



~ |2 



gl + gl 



0|2 



0|2 



(A33) 



where the dominant contribution Vq is given by Eq. (|9]). 



Appendix B: Bosonic Thermal Masses 



In order to calculate the daisy resummation Eq. (15) we require the thermal mass correc- 
tions lib. For the Higgs and singlet fields we compute the thermal mass corrections from the 
thermal effective potential using the procedure explained in this section. For the left-handed 
sneutrinos we use 
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which can be calculated by assuming that all species that are summed in AV^ are light. For 
the remaining bosonic species, we use the thermal mass functions calculated for the nMSSM 
by [i9\ in which the authors assumed that the Higgs, Higgsinos, electroweak gauginos, and 
SM particles were light. 

First, we evaluate the thermal effective potential correction Eq. ( [T3| as a function of the 
eigenvalues of the field-dependent mass matrices listed in Appendix lAl Let m^j be the j^^ 
eigenvalue of the z*^ mass matrix with has rii associated degrees of freedom. By writing the 



traces as a sum over eigenvalues and using that rii < for fermionic species, Eq. (13) can be 
written as 



2^ 



i bosonic 
i fermionic 



(B2) 



In the high-temperature limit, ml- ^ the bosonic and fermionic thermal functions can 
be expanded as 



JF{y) 



(B3) 
(B4) 



plus field independent terms. Second, we define the high-temperature thermal potential 
correction by imposing a sharp cutoff at m?- = 2T^ to obtain 



7", high 



48 ^ 



2m^^- i bosonic 
rij I < rfifj i fermionic 
mf. < 2T2 



(B5) 



Third, we extract the thermal mass corrections by differentiating with respect to the Higgs 
and singlet fields. 



rp2 



Q2 '^'^^ 



2^2 



(B6) 



=0,T=100 GeV 



where 0j G {H^-, v^}. The derivatives are evaluated at the origin in field space such that 



n,^. is accurate in the high-temperature vacuum. Because the derivative in Eq. (-B6) has 
only weak field dependence, we expect this expression for 11^^ to be accurate even for our 
multistep phase transitions in which the singlets have VEVs before the EWPT. The value of 
T used in Eq. [BQ] only affects the location of the cutoff in Eq. (55). We have chosen the 
temperature T = 100 GeV to be at the appropriate scale for our phase transitions and such 
that does not vary discontinuously in the region of parameter space with first order phase 
transitions. Using this procedure we obtain Hjjo 0.11 T'^.Huo 0.40 T'^^Iiy^. ^ 0.20 
over the region of parameter space with phase transitions. 



Appendix C: Analytic Derivation of Parameter Space Boundaries 



The boundaries in Fig. [5] can be understood analytically. In this section, we derive 
expressions for each of the boundaries and discuss the parametric dependence. 
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At the interface of regions I and II, the electroweak vacuum develops a tachyonic direction 
at tree level and det M^^^ = 0. Since M^j^^ is an 8 by 8 matrix, it would not be useful to 
write out its determinant. Instead, we observe that the tachyonic direction is directed along 
{H^/vi, HI/V2, Vtiv,. ^ 1, = 0}. 

At the boundary between region Ila and III, the minima at Xoi2 and yooo are degenerate 
at one-loop. Note that this degeneracy cannot occur at tree level. To see why, write 



K) (foi2) - K) (yl 



000 J 



(CI) 



with A^o" 



Vo (:?oi2) - 
symmetry ensures AVg" 



and AV^f 



Vo(^ooo) - V'o(yooo)- The tree level (Zg 



K) (a^ooo) 

= 0. Using the minimization equations Eq. (llOl) we can write 



AK,' 



{[sl + cos^ 2/3 + 6A^ sin^ 2/3] i;^ > 



(C2) 



where we have also neglected terms suppressed by fj>/f and At one-loop order we 
calculate the difference in the potential as 



(X012) - (yooo) = AVi" + AV,' 



(C3) 



where is the one-loop, zero temperature effective potential and AVj"'^ are defined anal- 
ogously as above. We expect AVi to be nonzero and sensitive to the radiative corrections 
because the (^3)^ symmetry is broken to Z3. The terms responsible for this symmetry 
breaking are the superpotential term W 3 Aifj'ifg'^f and corresponding A-term in the soft 
SUSY-breaking Lagrangian. We calculate 



eAn^AV^" = 6m^h log 



m 



ch 



>3/2^2 



+ 4 



4 1 



, 4 1 "^-ff2 4 1 

+ ^^^2log^^-2^"^± log 



ml 



g3/2^2 



m. 



it 



m 



+ m 



H2 



± 



m 



H2 



ch 



(C4) 



For the sake of discussion, we can approximate the logarithms in the second term as order 



one numbers and eliminate the soft masses using Eq. (10) to obtain 



QA-n'^AV^ 



6m^h log 



ch 



2m^h (t^^ + 4a CSC 2/3 - 4) + 2AX^v'^m^ 



(C5) 



Since we are simply trying to estimate the parametric dependence, we can approximate 
AV^ ~ AVq. By requiring that the minimum at X012 is not deeper than the EWSB vacuum, 
we obtain the bound AV^ + AVq > 0, which is saturated at the interface of regions Ila and 
III. This bound disfavors large rrich and large a because of the —m^^^a'^ term in AVi- 

At the boundary where regions III and lib meet, the EWSB is degenerate with the origin 
in field space at one-loop. Neglecting the radiative corrections we can approximate the 
splitting as (xq) - 1^° (yooo) ~ Vq (xq) - Vq (yooo) = AVJf with 



AV,^ 



I 2 2 



[{al + al) cos^ 2/3 + sin^ 2 
K sin 2/3 cr sin 2/3 



6A 



+ 



27A4 



+ 



ch 



27A3 



(C6) 
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To prevent the origin from becoming the global minimum we require AVq > which favors 
larger rrich and smaller a. 

At the boundary between regions lie and III, the one-loop potential has degenerate min- 
ima at yoi2 and yooo- We can compute the splitting AV^'^ = Vq (^012) — Vo ivooo) by neglecting 
the radiative corrections to find 

AVq^ = ^miy {a sin 2/3-2). (C7) 

The condition that the EWSB minimum at yooo is absolutely stable requires AV'o'^ > which 
imposes the lower bound cr > 2 esc 2/3 ^ 3 for tan /3 = 2.6. Figure IS] shows that the IIc-III 



boundary also depends on nich contrary to Eq. (C7), but this is a result of the radiative 
corrections. 



Appendix D: Selecting a CP-Even Vacuum 

In this appendix, we show formally how a superpotential contribution AW that breaks 
Z3 weakly can be constructed to make the CP conserving vacuum to have the lowest en- 
ergy perturbatively in the absence of any explicit CP violating parameters. Consider the 
superpotential 

W = Wo + AW (Dl) 

where AW represents a irrelevant operator perturbation to renormalizable Wq. We then 
have 



E 



dW 



dWo 



+ 23? 



/9Aiy\ fdWn 



(D2) 
(D3) 



where 



and 



AVi = 23? 



E 



dWn 



d(j)i 



/dAW\ fdWo 



(D4) 



(D5) 



to leading order in AW which we will call 0{6). As usual. Wo and AW are holomorphic 
polynomials in fields. Considering 0j — )■ 0* as a representation of Z2 which we will call 2, we 
have 

fdAW\ fdWo^ 



being a representation of 



V 9(t)i J 



2') = R. 



(D6) 
(D7) 
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If we assume all the coefficients of W and AW are real, we can write 

'fdAW\ fdWo 



23fJ 



R®Z2{R) (D8) 



where Z2(Z2(-R)) = R. Hence, we see that AVi is a singlet under Z2. Given that AVi is a 
polynomial in aj = 3fJ0j and bj = '^(pj and since under Z12 : {aj — )■ aj,bj — )■ —bj}, we must 



have 



where Sm represents a basis of Z2 singlet polynomial functions composed of bj and is a 
basis of polynomial functions composed of aj. Note that here c^^ = Hence, given that 

the part of the effective potential not associated with AW had a minimum at = v{s) where 
s e { — 1,0,1} parametrizes the Z3 fundamental representation elements and &j|i?(o) = is 
the singlet element, the energy shift due to AW to 0{6) is 

Ap{s) ^ Yl ^^Ms) = E E E ^ImPkiiv, COS ('^) })SU{v, Sin ('^'j }), (DIO) 

i i k m ^ ' ^ 

Note that Ap(l) = Ap(— 1). Hence, we only need to determine whether Ap(l) — Ap(0) > 
to see if CP singlet has the lowest energy. Since c^^ oc sgnAVT, we can simply flip 
the sign of Ap(l) — Ap(0) by flipping the sign of AW if the original choice of sign gives 
Ap(l) — Ap(0) < 0. Of course, all of this is under the assumption that the potential is not 
destabilized by the nonrenormalizable operators such that the smallness of the perturbation 
order 5 is meaningful. Stability is generic if the nonrenormalizable terms are dominated 
by the perturbations in the superpotential since the superpotential contribution is positive 
definite. 
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